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Any-order pion interferometry formulas for fixed pion multiplicity events and for 
mixed events are given. Multipion Bose-Einstein correlation effects on the two-pion 
interferometry and source distribution are studied. It is shown that generalized 
pion interferometry formula should depends on pion multiplicity distribution. Pion 
condensate is also discussed in the paper. 



1 Introduction 

Multipion Bosc-Einstcin(BE) correlations has now arous ed g reat interests among 
P hysicists0fflIlilE0li0-0li@01@0ili. Among others, 
Lamjand Loci were the first to suggest BASER concept in high energy physics. 
ZajctJ was the first to use Monte-Carlo methods to study multipion BE cor- 
relations effects on two-pion interferometry for fixed pion multiplicity event. 
Pratt □ first suggested that multipion BE symmetrization effects may lead to 
BE condensate in high energy heavy-ion collisions. Detail derivation of multi- 
pion BE correlation effects on single particle jspectrum and two-pion interfer- 
ometry was given by Chao, Gao and I in Ref.H. Recently Zimanyi and Csorgo 
suggested a new class density matrix and studyjrmltipion BE correlations and 
wavepacket effects on two-pion interferometry Qa. In Ref . EL2I, I have derived 
any- ord er pion interferometry formula for a special class density matrix. In 
Ref. liH, generalized any-order pion interferometry formulas are given which 
depends not only on the correlator as assumed previously but also on the pion 
multiplicity distribution which was neglected in the previous studies. Multi- 
pion BE correlation effects on the source distribution was studied in Ref. Ell. 
This paper is arranged as follows: In section 2, any-order pion interferometry 
formula for fixed pion multiplicity event are given. In section 3, any-order pion 
interferometry formula for mixed events are derived. In section 4, methods on 
multipion BE correlation simulation are discussed. Conclusions are given in 
section five. 



2 Any-order pion interferometry formulas for fixed pion multiplic- 
ity events 

Assume the source is totally chaotic. It is easily checked that the n-pion 
inclusive distribution in n-pion event can be expressed as El 

n 

^(Pl.--'.Pn)=X)II^>W) (1) 

a j = l 



Pij = I d 4 xg(x, Pl+P: > )exp(i(pi - p 3 ) ■ x). (2) 



with 



Here g(x, k) is a Wigner function which can be explained as the probability of 
hiding a pion at point x with momentum p. cr(J) denotes the j-th element of 
a permutations of the sequence {1,2,- • •, n}, and the sum over a denotes the 
sum over all n! permutations of this sequence. Then the normalized k pion 
inclusive distribution in n pion events can be expressed as 
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It is easily checked that P„(pi, • • •, pk) can be re- written as 
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with 



i=l 

and 



1 - /" 
u(n) = -^2«i(n-i)J d P Gi(p,p), w(0) = l (5) 



Gi(p,q) = y /o(p,pi)rfpip(pi,p 2 ) • • • dpi_ 1 p(p i _ 1 ,q). (6) 

Here er(i) denotes the z-th element of a permutations of the sequence {1, 2, • • 
•,/c}, and the sum over a denotes the sum over all fc! permutations of this 
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Figure 1 : Multipion BE correlation 
effects on two-pion interferometry. 
R=3.5fmandA=0.14GeV. The 
total momentum of the two pions 
is zero 
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Figure 2: Multipion BE correlation 
effects on source distribution in 
momentum space. R=3.5fm and 
A = 0.16 GeV. 



sequence. Then the fc-pion (k < n) interferometry formula for n-pion evnets 
can be expressed as 



C£(pi,- • -,p k ) 



^(Pir -,Pk) 

n-=i p «(pj) ' 



We assume the source distribution as 



g(r,t,p) = n • ( 



f /2 exp{- 



2A 2 



(7) 



)• (8) 



Here no is a parameter. Then the multipion BE correlation effects on two-pion 
interferometry for fixed evnets are shown in Fig.l. It is easy to see that as 
pion multiplicity increases, multipion BE correlations effects on the two-pion 
interferometry effects becomes stronger. In the limit n^oowe have C\ = 1. 



3 Any-order pion interferometry formula fo mixed events 

In the experimental analyses, one normally mixed all events to study pion 
interferometry. The fc-pion inclusive distribution can be expressed as 

oo 

Wfc(pi, ■•■ lPk ) = p B -n-(n-l)-(n-H l)^„ fe (Pi, ■ ■ ,Pk). (9) 

n—k 

Here P n is the normalized pion multiplicity distribution. It is easily checked 
that 

f N k ( Pl ,--;p k )f[d Pj = (n(n-l)---(n-k + l)) , (10) 
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then the /c-pion interferometry formulas can be expressed as 
n , \ N k (px,-- -,p k ) 

n j= i^i(Pj) 

The discussion about the definition of two-pion interferometry formulas can 

lool 

be found in the paper by Miskowiec and Voloshinl££l and can also be found in 
Ref.El. Using Eq.(4) , the /c-pion inclusive distribution can be expressed asliil 
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It is interesting to notice that if P n = w(n)/ J) as assumed in 

ve have the following modified k pion inclusive distribution^: 

N k (pi, ■ ■ -Pk) = E ^Mi) ' ' ' H ^(k) (I 4 ) 



with 

oo 

%=lf(p i ,Pj) = ^G n (p i ,p j ). (15) 

n=l 

One interesting property about Eq.(15) is: Eq.(15) is very similar to the pure 
n pion inclusive distribution (Eq.(l)) but with a modified source distribution 
S(x, K) which satisfies 

ff(Pi,Pj) - j S(x,^±^)exp(t( Pl -p j )-x)d i x. (16) 

Based on Eq.(16), we will study multipion BE correlations effects on the 
source distributionsEll. From Eq.(16), the single particle spectrum distribution 
,Pi™^(p), can be expressed as: 

W (p) = = jS{x p)*x = j = j s(X)p)dixdp _ 

(n) (n) J J 

(17) 




Similarly the source distribution in coordinate space ,Pi (r) , can be expressed 
as 

^ (r) = j5fep)dpdt [ (18) 

Assume the input g(x,p) as Eq.(8), multipion BE correlation effects on the 
source distribution are shown in Fig. 2 and Fig. 3. It is clear that as (n) in- 
creases, multipion BE correlations effects on the source distribution becomes 
stronger. Multipion BE correlation make pions concentrate in momentum and 
coordinate space. It is easily checked that the half width of the source distri- 
bution in momentum space and coordinate space A e /j and R e ff satisfies the 
following relationship: 

2 ^<A e// <A, ^<R eff <R. (19) 

If we take ^g^s — > oo. we have Reff = \J 53 an d &eff — \l 55- The 
source size satisfy the minimal Heisenberg relationship. That is all pions are 
concentrated in a single phase space and pion condensate occurs. But the pion 
multiplicity distribution is Boson form:P„ = (^"^n.+i ■ 

In general, pion interferometry formula depends not only on the correlator 
Pi j but also on P n as shown in Eq.(12). The later property was neglected in 
the previous studies. Recently, Egger, Lipa and BuschbeckH results seems 
support the conclusion presented here. Previous studies based on pure multip- 
ion interferometry formulas have shown that the difference b etwe en three-pion 
interferometry and two-pion interferometry are not so larger c^E^I. Certainly, 
there is also another possibility that the there are some correlation among the 
emitted pions which was not include in the above derivation. So one needs 
more studies about the two-pion and higher-order pion interferometry! 

If we assume the pion state as 



^a„|n>. (20) 



Here a n is a parameter which connected with pion multiplicity distribution. 
\n > is the non- normalized pion n pion state which can be expressed as 

( /rf p j( p) a+ (p))" 
n\ 

Here a + (p) is the pion creation operator and j(p) is the pion probability 
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amplitude. Then the pion multiplicity distribution P n can be expressed as 

Pn = ' a "' 2 ; ( " ) y (22) 

The two-pion interferometry results for different P n (different a„)is shown in 
Fig. 4. One can see clearly that there are big differences among the two-pion 
interferometry results for different P n . 



4 Monte-Carlo simulation of multipion interferometry 

To put two-pion BE correlations in event generator has been studied by dif- 
ferent authors c3. To put multipion BE correlations in event generator is a 
very difficult task as shown in Ref. Recently Fialkowski et alLLil used 

the methods of Bialas, KrzywickiliB and Wosiekliil and implemented approxi- 
mately multipion BE correlation in JETSET/PYTHIA to study multipion BE 
correlation effects on the W mass. Due to the fact that there is no space-time 
picture in the model JETSET, So they had to put pij by hand. On the other 
hand for model which has space-time picture, one can construct g(x,p) accord- 
ing to following method H'Q: D(x,p) is a total classical function which is the 
output of the evnet generator and can be written down as: 

D(x,p) = Y / 6(x-x(i))6( P -p(i)). (23) 



Here x(i) and p(i) are the i particle coordinate and momentum. The natural 
and easiest way to construct g(x,p) the semiclassical function is to replace the 
above delta function with a Gaussian wavepacket 

*(y, fc)=EII exp(- fa-f^ 2 ) exp(-^|% (24) 

In general Sxj {Spj) should not be the same. But for simplicity one can set 
them the same value and take Spj = The same discussion can be found 
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in the Ref. Ifrl. Unfortunately the calculating work for the above methods to 



put multipion BE correlation in event generator is so large that a new method 
which enable us to implement quickly multipion correlation in event generator 
is a debate for physics in this field. 
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Figure 3: Multipion BE correlation 
effects on the source distribution 
in coordinate space. R=3.5 fm 
and A=0.16 GeV. 
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Figure 4: Two-pion BE correlation 
for different pion multiplicity distri- 
bution. R=3.5 fm and A=0.1 6GeV. 
The mean pion multiplicity for 
above three cases are 1 7.8. 



5 Conclusions 

In the above, we mention nothing about the effects of flowll^ and resonancelll 
on the pion interferometry. As our main interests is the BE statistical effects 
and the flow and resonance can be included in the function g(x 7 p). Certainly 
the multipion coulomb effects is still a open question in this field. Energy 
constraint effects on the pion interferometry can be studied by the method 
presented in Ref.LJ. In this paper, we have derived any-order pion interferom- 
etry formula for fixed multipion events and for mixed events. In general pion 
interferometry formulas should depends on not only on the correlator but also 
on the pion multiplicity distribution which was neglected in previous studies. It 
is shown that as the pion multiplicity becomes very larger all pions concentrate 
in a single phase space and pion condensate occurs. But the pion multiplicity 
distribution is BE form. The methods to put multipion BE correlations in the 
event generator is discussed. 
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